Effects of non-standard interactions in the MINOS experiment 

Mattias Blennow|*] Tommy Ohlssonj^ and Julian Skrotzkjf] 
Department of Theoretical Physics, School of Engineering Sciences, 
Royal Institute of Technology (KTH) - AlbaNova University Center, 
Roslagstullshacken 21, 106 91 Stockholm, Sweden 

00 ■ We investigate the effects of non-standard interactions on the determination of the 

o ; 

. neutrino oscillation parameters Am|^, 023> and ^13 in the MINOS experiment. We 
CN ■ 

^-^ ' show that adding non-standard interactions to the analysis leads to an extension of 

! 

^ ' the allowed parameter space to larger values of Amg^^ and smaller ^23) and basically 

, removes all predictability for ^13 . In addition, we discuss the sensitivities to the non- 

\ standard interaction parameters of the MINOS experiment alone. In particular, we 

' examine the degeneracy between ^13 and the 

^ ■ 

■ We find that this degeneracy is responsible for the removal of the ^13 predictability 

(N : 

, and that the possible bound on |eer I is competitive with direct bounds only if a more 

■ 

I stringent external bound on ^13 is applied. 

ix. 

^ ■ I. INTRODUCTION 

> : 
•1—1 , 

^ \ In the years that have passed since neutrino oscillations were first observed at the Super- 



Kamiokande experiment in 1998 [1], there has been a remarkable experimental development 
in neutrino physics. For example, from the results of solar 



m 



3|, IJ] and long-baseline reactor 
neutrino experiments, we now know that the solution to the solar neutrino problem is 
given by the large mixing angle (LMA) solution with Am^i — 8 • 10~^ eV^ and 612 — 33.2°, 
and from atmospheric jl| and accelerator 3, Isl, Q] neutrino experiments, we know that 
lAmg^l ^ 2.5 ■ 10^'^ eV^ and that 6'23 is close to maximal (i.e., 6'23 = 7i"/4). In addition, 
analyses of the L/ E binned Super-Kamiokande and KamLAND [6.] data even show the 
oscillatory behavior of the neutrino flavor conversion probability. 
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The fact that neutrino oscillations occur implies that neutrinos have non-zero masses, 
which requires physics beyond the standard model of particle physics (SM). Thus, neutrino 
physics seems to be a viable window to explore physics beyond the SM. A feature of m any 
extensions of the SM is the existence of non-standard interactions (NSI) (see, e.g., Ref. 11| 
for a recent review) between neutrinos and other fermions, including the first generation 
fermions which make up most of the matter that we experience in everyday life. In particular, 
effective four-fermion operators arising from such NSI will inevitably affect the dispersion 
relations for neutrinos propagating in matter through coherent forward scattering similar 
to that of the Mikheyev-Smirnov-Wolfenstein (MSW) effect 12, Q, ll^, which is usually 
considered in neutrino oscillation analyses and responsible for the conversion of solar z/g into 
Ufj^ and u-r. With new generations of neutrino oscillation experiments in the planning stages, 
we expect to probe the yet unknown parts of the parameter space for neutrino oscillations 
and to decrease the experimental uncertainty in the parts where we have only pinpointed 
certain regions. In such precision experiments, it may happen that even small contributions 
of NSI to the matter effects can play a role in distorting the measurements of the standard 
neutrino oscillation parameters or, more excitingly, that NSI can even be observed through 

nn 

the very same effects [15|,116|. 

The Main Injector Neutrino Oscillation Search (MINOS) jol is an accelerator based neu- 
trino oscillation experiment with a baseline of 750 km reaching from Fermilab, Illinois to 
the Soudan mine, Minnesota in the United States. It is an experiment designed to mea- 
sure the neutrino oscillation parameters Am|^ and 623, but it may also improve the bound 
on the leptonic mixing angle 613. In this paper, we discuss the implications of including 
NSI in the analysis of the MINOS data. We focus on how the introduction of NSI affects 
the experimental bounds on the standard neutrino oscillation parameters, but also discuss 
what bounds MINOS itself could put on the parameters Sa/s, which describe the NSI on a 
phenomenological level. 

Non-standard interactions in the MINOS experiment have been previously studied by 
Friedland and Lunardini in Ref. While they focus on constraints which are put by 

the combination of MINOS and atmospheric neutrino oscillation experiments, we focus on 
the constraints that can be inferred from the MINOS experiment alone. In particular, we 
consider the z/g appearance channel and its implications for the leptonic mixing angle 

6'i3 and the effective NSI parameter Ser in detail. Also in Ref. [l8|], the effects of NSI on 
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the z/g appearance channel at MINOS were studied, focusing on the oscillation probability 
P^e- One of the conclusions of Ref. 3] was that, in the most optimistic case, the oscillation 
probability will be so large that it cannot be described by the standard neutrino oscillation 
scenario alone, and thus, implying the existence of NSI. Our numerical simulations will 
show that l^erl — 2.5 (which is above the current experimental bound) would be needed to 
establish NSI unless further external constraints can be put on sin^(20i3). 

The paper is organized as follows. In Sec. [Tll we introduce the framework of neutrino 
oscillations including the effects of NSI. Section UTTl deals with analytic considerations for the 
neutrino oscillation channels relevant to the MINOS experiment, while Sec. IIVI presents the 

n n 

results of our numerical treatment using the GLoBES software (JJ, l20(. Finally, in Sec. |Vl 
we summarize our results and give our conclusions. 



II. NEUTRINO OSCILLATIONS AND NSI 



In this paper, we will use the standard three-flavor neutrino oscillation framework with 
an effective vacuum Hamiltonian given by 



Ho = — f/diag(0, Am2„ Am2j[/t 



(1) 



in flavor basis. Here E is the neutrino energy, AmJ^_= mf — m] are the neutrino mass 



squared differences, U is the leptonic mixing matrix 21 1 
/ 

U = 



C13C12 



C13S12 



S12C23 - Ci2S23Sl3e'^ C12C23 - Si2S23Sl3e'^ S23C13 



v 



'Sl2'523 — Ci2C23Si3e 



i<5 



"Cl2"S23 ~ Si2C23Si3e''^ C23C13 



(2) 



Cij = cos{9ij), Sij = sm{9ij), 9ij are the leptonic mixing angles, and 6 is the CP- violating 
Dirac phase. In addition, the standard matter effect on neutrino oscillations is implemented 
through the effective contribution 12, 1^ 



i^MSW = diag{V2GFNe, 0, 0) = diag(l, 0, 0) 



(3) 



to the vacuum Hamiltonian, where Gp is the Fermi constant and N^. is the electron number 
density. 

We are interested in examining the effects of introducing NSI between neutrinos and other 
fermions that reduce to effective four-fermion interactions. These NSI can be described by 
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a Lagrangian density of the form 

G 



C 



NSI 



^ E ^5fe^(l-75)^/3][/7.(l + «75)/], 



(4) 



a=±l 



where the e(^^ give the strength of the NSI. In analogy with the MSW effect, terms of this 
type will give an effective contribution to the neutrino oscillation Hamiltonian, which will 
be of the form 



NSI 



V 



^ee ^er 



\ 



F F F 
^e^( '^ii.ji. ^f^T 



c c c 
y ^er ^fiT ^rr J 



(5) 



where 



Nf is the number density of fermions of type /, and we have assumed an unpolarized 
medium. For bounds on the parameters e-^^, see Refs. 22|, [23]. Generally, the NSI involving 



are quite well constrained, while the bounds on the other NSI {i.e., Fee, £eT, and Frr) 
are of order unity. Thus, we will focus on NSI which do not involve z/^ interactions. In 
the remainder of this paper, we will work with the effective parameters Fai3, assuming them 
to be constant, which is a good approximation as long as the matter composition does not 
change significantly along the neutrino baseline. The full Hamiltonian is then given by 



H = Ho + H, 



MSW 



(6) 



Effects of this type have beeri previous 
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y studied in Refs. [12, 
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III. ANALYTIC CONSIDERATIONS 



In this section, we present some analytic considerations which are valid mainly for weak 
NSI. These will prove useful in understanding the numeric results in the next section. 

The main objective of the MINOS experiment is to measure the neutrino oscillation 
parameters Am|^ and 6*23. The neutrino oscillation channel used is the z/^ disappearance 
channel, which is sensitive to the survival probability P^^. The leading terms in the 
expression for P^^ are 



P„„ ~ 1 - sin2f2^o,) sin^ ( —J^L 



4E 



(7) 
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where three-flavor effects due to Aml^ and ^13 have been neglected. Equation ([7j) can be 
easily derived using the effective two- flavor Hamiltonian of the v^-Vr sector, i.e., 

^ Amfi / - cos(2^23) sin(2^23) \ 
° \^ sin(2^23) cos(2023) / ' 

For the base-line and matter potential relevant to the MINOS experiment, the off-diagonal 
NSI parameters will not be sufficient to introduce large transitions to Ve-, and therefore, the 
NSI can be effectively discussed in the same two-ffavor framework, where their contribution 
to the effective neutrino oscillation Hamiltonian is 

<si = Ve^^ [ I 1 ] ■ 

The stringent way of treating this situation is to introduce the parameter E^f clS Qb perturba- 
tion. It is then easy to show that the two-ffavor approximation holds for l^erP^^-^^ ^ or 
l^g^p ^ 5.8 in the case of the MINOS experiment. We note that neutrino oscillations with 
the effective Hamiltonian H^^ = Hq^ + if^si equivalent to standard two-ffavor neutrino 



oscillations in matter with the substitutions Am^ ^^31' ^ ^23) and V 



which we know that the effective neutrino oscillation parameters are given by 



12 



13 



or 



Am' = AmU, sin^(2g) = '^^'^f , (10) 



where 

2EV 

+ C0S(26'23) 



+ sin2(2^23). (11) 



Thus, for a fixed energy E, it is always possible to choose Err in such a way that C, = sin(26'23), 
leading to Am^ = sin(2^23)Am|| and sin^(2^) = 1. Therefore, we expect, when including 
the effects of NSI, a degeneracy between the standard oscillation parameters Am|]^ and 
sin^(2^23) and the NSI parameter Err, i-^-, if we measure Am^ = AttIq and sin^(2^) = 1, 
then this could just as well be produced from a smaller mixing angle and larger mass squared 
difference by the effects of the NSI. The fact that different Err will be needed in order to 
reproduce this effect at different energies implies that this degeneracy can be somewhat 
resolved by studying the neutrino oscillation probability at different energies (as in the case 
of an actual neutrino oscillation experiment measuring the neutrino energy, e.g., the MINOS 
experiment). However, if the energy range is not broad enough, then the degeneracy will 
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sin (2e„) = 1 



sin^ie^j) = 0.75 

5111^2623) = 0.5 



FIG. 1: The analytic result for the neutrino oscillation probability P^^ as a function of the neutrino 
energy E for different values of sin (26I23) and Am|i on the NSI de generacy involving (only 
the value of sin^(2^23) is displayed in the figure) in the two- flavor scenario. The NSI parameter 
Err has been chosen in order for the effective parameters to coincide aX E = 2> GeV. 



still manifest itself in the form of an extension of the sensitivity contours when including NSI 
into the analysis. In Fig. [H we show the neutrino survival probability P^^ as a function of 
energy E for different choices of sin^(26'23) and Amg^ which are on the NSI degeneracy. Note 
that this figure is only provided for illustrative purposes in order to show the degeneracy 
and that the values of Srr needed for sin^(2^23) = 0.5 and 0.75 are relatively large (about —7 
for sin^(26'23) = 0.5). However, the degeneracy does not need to be exact in order to extend 
the sensitivity countours and also smaller values of Srr will be enough for this purpose. 



In Ref. 



161], it was shown that, including the first-order correction in ^q,^, the ef- 



fective three-flavor mixing matrix element t/es is given by (to zeroth order in the ratio 

Ue3 = Ue3 + ^er-Z 2~^23! (12) 

which is valid as long as the individual contributions remain small (from the CHOOZ bound 
42] Ue3 is known to be small and the absolute value of the NSI contribution is of the order 
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of 0.2 for l^eT-l = 1). As the oscillation probability P^e is expected to be sensitive to the 
effective mixing angle ^13, this will lead to an additional degeneracy between the parameters 
9i3 and Eer- The effects of the other NSI parameters are suppressed by a or S13. 



IV. NUMERIC SIMULATIONS 

For our numeric simulations, we used the GLoBES software [lol, [2^ which was extended in 
order to accommodate the inclusion of NSI. The Abstract Experiment Definition Language 
(AEDL) files used to describe the MINOS experiment were modified versions of the MINOS 
AEDL files provided in the GLoBES distribution and they were based on Refs. 43|, |4J, 
451]. These AEDL files correspond to a MINOS running time of five years with 3.7 ■ 10^° 
protons on target per year. The neutral- and charged-current cross-sections were taken from 
Refs 



46 



471] as provided by the GLoBES distribution. 



The disappearance and appearance channels were simulated in a neutrino energy interval 
of 1-6 GeV, since the majority of the neutrinos in the NuMI beam are in this range. For the 
simulations, the neutrino energy interval was binned into 30 equal bins. The matter density 
was assumed to be constant with a value corresponding to the matter density of the Earth's 
crust, i.e., V = 1/1900 km~^. The simulated neutrino oscillation parameters are shown in 
Tab. [11 The choices for Amg^ and sin^(2^23) are inspired by the preliminary MINOS results 
Q], which are almost equivalent to the K2K results Q], and the simulated values of all NSI 
parameters are zero, in order to possibly obtain useful sensitivities for NSI detection. In 
all simulations, we have used the full numeric three-flavor framework and the parameters 
not presented in the figures (including NSI parameters such as the phase of egr) have been 
marginalized over unless stated otherwise. It should also be noted that normal mass hierar- 
chy, i.e., Amg]^ > 0, was assumed for the simulations shown. The results for inverted mass 
hierarchy are similar and no distinction can be made between the hierarchies. Furthermore, 
sin^(26'i2) and Am^i were kept fixed at the values given in Tab. [T]for the simulations, since 
MINOS is not sensitive to these parameters. In addition, the CP-violating phase 6 was kept 
fixed for the simulations of the parameters governing the disappearance channel, since its 
effect in this channel is small. However, for the simulations of the parameters governing 
the appearance channel {i.e., 613 and Eer), we marginalize over S, since it is important to 
include the effects of a possible relative phase between Ues and Est- The explicit choice of 
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sin2(26'i2) = 0.8 


Amj^ = 7 • 10-5 eV2 


sin2(26'i3) = 0.07 


Am|i = 2.74 • 10-3 eV^ 


Sin2(2023) = 1 


6 = - 



TABLE I: The neutrino oscillation parameters used in the simulations. 

7r/2 for the simulated value of 5 does not affect the results of our simulations significantly. 
In all figures, we show the combined results of the disappearance and appearance channels. 
Furthermore, when the standard neutrino oscillation parameters are marginalized, then we 



assume 20 % external error (la) for Arrig^ and 6*23 and an external error of O.OGvr for 9 
(corresponling roughly to the CHOOZ bound). 



nal errors in accordance with direct bounds 
atmospheric neutrino oscillations 



22 



or the NSI parameters, we assume exter- 



40 



23( 1 and with the results of high-energy 



4l| . The high-energy sample of atmospheric events 
indicate that muon neutrinos oscillate also at higher energies. If the second eigenvalue A2 
of the matter interaction part of the full Hamiltonian (including both NSI and the standard 
matter effect) is too large (such that IA2I ^ Am\^/ {2E)), then matter effects will entirely 
dominate the neutrino flavor propagation. Since we have assumed e^a = ^p^i = 0, this would 
mean that muon neutrinos would be fully decoupled in contrast to experiments. The result- 
ing constraint in the NSI parameter space has the shape of a parabola in the eT-^-l^erl-plane 



as A2 = corresponds to Ej 



£ee) |4Q|, |4l|. However, the study of high-energy 



events of different flavors (see, e.g., Ref. 48] ) may provide further information such as the 
composition of the state which muon neutrinos oscillate into, which in turn will be related 
to the NSI parameters. The atmospheric constraints were implemented by setting a prior 
of IA2I < 0.2V for the second eigenvalue. The results are not particularly sensitive to the 
specific prior chosen. 



A. Degeneracy of Stt 

Figure [2] shows the predicted sensitivity limits of the MINOS experiment (according to the 
experimental setup given above) in the sin^(26'23)-Am3]^ plane with and without the inclusion 
of NSI. From this figure, we can clearly observe the extension of the sensitivity contours 
according to the discussion in the previous section. The reason why the contours do not 
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99 % C.L. 
95 % C.L. 
90 % C.L. 
99 % C.L. 
95 % C.L. 



(no NSI) 
(no NSI) 



90 % C.L. (no NSI) 
NSI degeneracy 
Best-fit 




0.6 



0.7 



0.8 



0.9 



sin (2023) 



FIG. 2: The sensitivity limits in the svc? {292'i)-~^m^^ plane (2 d.o.f.) for the combined appearance 
and disappearance channels. The colored regions correspond to the sensitivities when including 
NSI, while the black curves correspond to the sensitivities under the assumption that NSI are 
negligible. The dashed blue curve marks the NSI degeneracy involving Err where Amg^^ sin(2023) = 
2.74 • 10"'^ eV'^. The best-fit point corresponds to the parameter values used in the simulation. 



extend to sin^(2^23) = is based on the fact that the MINOS experiment is not using a single 
neutrino energy, but rather has a continuous energy spectrum. For a fixed ^ = sin(2^23) 
will only be fulfilled for one specific energy and the effective neutrino oscillation parameters 
will become energy dependent. Although ^ = sin(2^23) niay still be approximately fulfilled 
in some finite energy range, for lower sin^(2^23), the energy dependence will become strong 
enough for the MINOS experiment to detect it and this is where the sensitivity contours in 
Fig. [2] are cut off (c./.. Fig. [Tj). In addition, improved external bounds on the NSI parameter 
Err could lead to a cutoff for the extended sensitivity contours. 



B. Degeneracy of e^T 

In Fig. [3], we show the predicted sensitivity in the sin^(26'i3)-|£:eT| plane for one degree 
of freedom (signifying that we consider the sensitivities for the two parameters separately). 
This figure has been constructed assuming NSI parameters along the parabola allowed by 
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- 95 % C.L. 

90% C.L. 

4 Best-fit point 
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0.1 

sin\2e|3) 



0.15 
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sin'(2e„) 



FIG. 3: The sensitivity limits in the sin^(20i3)-|eeT| plane (1 d.o.f.) for the combined appearance 



and disappearance channels. The left panel corresponds to a simulated sin^(20 



13 J 



and the 



right panel to sin^(20i3) = 0.08. The black curves correspond to |?7e3p [c/-, Eq. (fT2|) ] equal to the 



simulated value of sf^ for E = 2.3 GeV. In this figure, See = and is chosen along the parabola 



allowed by atmospheric neutrino experiments 



40 



4li. 



atmospheric neutrino experiments 



40 



4l| and £ee = (allowing for general values of Eee 
slightly extend the contours). The results are shown for both sin (2^13) = and sin^ (2^13) = 
0.08. As can be seen from the left panel of this figure, the MINOS experiment is sensitive 
to sin^(2^^i3) which is a factor of two below the CHOOZ bound if we do not take NSI into 
account. However, if we include the effects of NSI, then the bound put on sin^(20i3) will 
depend directly on the external bound on l^erl- Already for a bound of |eer| ^ 0.5, the 
bound that MINOS is able to put on sin^(2^i3) has deteriorated to the CHOOZ bound and 
is quickly getting worse for less stringent limits on {eerl- We can also clearly observe the 
sin^(26'i3)-|£er| degeneracy discussed in the previous section. The sensitivity contours of 
Fig. [3] contain the degeneracy curves corresponding to Ue3 equal to the simulated values of 
5^3 for neutrino energies in the MINOS energy range. For possible bounds on l^erl, we need 
to consider some external bound on sin^(26'i3). With the current CHOOZ bound, the bound 
that could be put by the MINOS experiment is l^erl ^ 2.5, which is not competitive with the 
current direct bounds on the specific NSI. If the sin^(2^i3) bound is improved by an order of 
magnitude {e.g., by future reactor experiments 



49 



50 



51 



52 



53|), then the MINOS bound 



on l^erl could be improved to Iserl ^ 1, which is still of the same order of magnitude as the 
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direct NSI bounds. Thus, it seems that in order to put constraints on NSI from neutrino 
oscillation experiments, we would need an experiment with better sensitivity than MINOS. 

In the right panel of Fig.[3l we show the results for sin^(2^i3) = 0.08 (s^g = |f/e3p — 0.02). 
Again, we can observe the degeneracy in the sin^(2^i3)-|eer| plane along the direction of 
|f/e3p ~ 0.02. In this case, MINOS will be able to tell us that f/es is non-zero, but not 
whether this is the effect of non-zero ^13 or non-zero ^er (or a combination). Because of 
this degeneracy, the result is very similar if considering no n- zero e^t while 6*13 = 0. When 
considering and Ser which are simultaneously non-zero, the resulting sensitivity contours 
depend on the magnitude of the simulated f/es [e.^^., if the two terms in Eq. (fT2|) cancel, then 
we obtain sensitivity contours similar to the left panel]. 



C. Prospects of detecting NSI at MINOS 

In Fig. m we show the prospects of detecting NSI at the MINOS experiment, i.e., we 
show the regions of the NSI parameter space where the sensitivity contours do not contain 
the standard oscillation scenario {Sap = 0). As can be seen from this figure, the current 
situation is such that MINOS will not be sensitive to any of the NSI parameter values which 
are not already excluded (the excluded regions are Eee > 2.6, Eee < —4, and l^erl > 1-9 



in accordance with Ref. [22]). Increasing the MINOS running time to 15 years does not 
improve significantly upon this result. However, if the external bounds on the standard 
neutrino oscillation parameters are improved by a factor of four, then MINOS will be able 
to detect Iserl of about one at a confidence level of 90 %. This is due to the improved 
bound on ^13 leading to a breaking of the sin^(2^i3)-|eer| degeneracy. The improvement in 
the external bound on ^13 used in the figure corresponds to an upper limit of sin^(2^^i3) ~ 
10~^, which is slightly below the Double Chooz sensitivity limit [54] (note that the reactor 
neutrino experiments are not sensitive to NSI, since they operate at very low energies). The 
asymmetry of the figure with respect to Err = is a result of how Err affects the effective 
[/e3. Since we are not in the perturbative regime, we treat the case with ^13 = exactly 
in Err and perturbatively in E^r- The result of this is an equation similar to Eq. f|T2|) but 
where 623 is an effective quantity which depends on Err- The main influence on the NSI 
sensitivity is given by the appearance channel, with no or only a small contribution from 
the disappearance channel. 
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FIG. 4: The prospects for detecting NSI at the MINOS experiment. The sohd curves correspond 
to our default setup with five years of running time and external bounds on Am^^, 623, and ^13 
as defined in the beginning of this section. The dashed curves correspond to an increased running 
time of 15 years, whereas the dotted curves correspond to five years of running time, but with 
external bounds which have been improved by a factor of four. The NSI parameter See has been 



With the advent of new precision measurements of the neutrino oscillation parameters, 
it is important that we understand the phenomenology of the physics that could affect these 
measurements and give rise to erroneous interpretations if not taken properly into account. 
In addition, putting constraints on such physics from neutrino oscillation experiments alone 
is also an intriguing idea. In this paper, we have studied the influence of including NSI into 
the analysis of the MINOS experiment by analytic arguments and by using the GLoBES 
software in order to simulate how the sensitivity to the ordinary neutrino oscillation para- 
meters is affected by the introduction of NSI. We have also studied the prospects of putting 



chosen along the parabola allowed by atmospheric neutrino experiments [40|, |41|. 




V. SUMMARY AND CONCLUSIONS 
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bounds on the effective NSI parameters directly from the MINOS data. 

Our analytic results show that the disappearance channel (z/^ —>■ oscillations) is mainly 
affected by the effective NSI parameter Erri while the appearance channel (z/^ Ve) is 
mainly sensitive to Ser- The effect of including NSI into the analysis of the disappearance 
channel is that the sensitivity contours are extended to larger Am|^ and lower sin^(2^^23) by 
the introduction of a degeneracy due to Ett as described in the analytic treatment. These 
analytic considerations are supported by our numerical simulations, which are performed 
using the full three-flavor framework. 

In the numeric analysis of the appearance channel, the degeneracy b 
mixing angle 6'i3 and the effective NSI parameter Eer described in Refs. 
difficulties in placing bounds on either of these parameters unless a stringent bound for the 
other parameter is imposed by external measurements. With an external bound on \EeT\ of 
the order of 10~^, the MINOS experiment would be sensitive to values of sin^(26'i3) down to 
about 0.07, to be compared with the present bound of approximately 0.13 from the CHOOZ 
experiment. However, this sensitivity rapidly deteriorates with less stringent bounds on 
|eer| and the sin^(20i3) sensitivity for the NSI parameters, which are phenomenologically 
viable today, is clearly worse than the CHOOZ bound. On the other hand, if external 
measurements show that sin^(20i3) < 0.01, then the MINOS experiment should be able to 
place a bound on \E(.t\ of the order of unity, which is of the same order as the present bounds 
from interaction experiments. With the current CHOOZ bound, the bound that could be 
put on l^erl from the MINOS experiment is about a factor of 2.5 larger than this, even with 
an increased running time of 15 years. However, a signal in the MINOS appearance channel 
would indicate that either ^13, Eer-, or both are non-zero as this would imply U^z 7^ 0. 

It should also be noted that the results presented here do not depend on the neutrino 
mass hierarchy {i.e., if Amg^^ > or Amg^^ < 0), or whether or not we try to make a fit with 
the same mass hierarchy as the one used in the simulation. 

In conclusion, it seems that the MINOS experiment is very close to being able to put 
a useful bound on \Eer\ if sin^(2^i3) could be further constrained by, e.g., future reactor 
experiments. Thus, the next generation of neutrino oscillation experiments should be able 
to put bounds on |ee^| which are more stringent than the direct bounds from interaction 
experiments. 
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